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Abstract All engineers who work to improve the work efficiency of the a Stirling 
engine try to replace the crank mechanism of compression cylinder by another one. 
It could be useful a mechanism able to perform a dwell during the whole period of 
expansion. This stop of the compression piston has to enlarge the active area of 
the PV diagram. As consequence it is expected an efficiency improvement. This 
request could be very well accomplish by a cam mechanism. But this one is too 
large, too complicate, expansive and unreliable for solar systems as shown in [1]. 
Hence, the paper tries to find out a linkage able to perform an equivalent dwell as a 
cam mechanism. The purpose of the analysis is to find out the movement functions 
yC(t) and yE(t) to reveal if there are in a proper phase relation between them or not. 
The proper phase relation demonstrates that the mechanism is good for the applica¬ 
tion. Second, it must show that yE(t) has a large enough dwell, in accordance with 
the expansion stroke time of the expansion cylinder. 


Keywords Bar mechanism • Crank rod mechanism • Dwell mechanism • 
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1 On the Working Cycle of Stirling Engine 

The PV diagram of a-Stirling engine it is far away from an ideal PV diagram. The 
principle schema of a a-Stirling engine is shown in Fig. 1. As seen the engine has 
two crank rod mechanisms, one for the compression cylinder and the other for the 
expansion cylinder. The movement functions for both pistons are sinusoidal, as 
shown in Fig. 2. 
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Fig. 1 The standard a- 
Stirling engine 




crank angle 


Fig. 2 Piston movement functions 

Therefore, all the engineers who struggle to improve the work efficiency of the 
engine try to replace the crank mechanism of the compression cylinder by another 
one. It could be useful a mechanism able to perform a stop during the whole period 
of expansion. This stop of the compression piston has the role to enlarge the active 
area of the PV diagram shown in [2] and as a consequence it is expected an effici¬ 
ency improvement in [4]. This request could be very well accomplish by a cam 
mechanism. But this one is too large, too complicate, expensive and unreliable. 
Hence, we try to find out a bar mechanism able to perform an equivalent stop as a 
cam mechanism. 
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In order to reach the dwell displacement function, without using a cam mecha¬ 
nism, the authors try to create a new bar mechanism. 


2 A New Dwell Mechanism Dedicated to a-Stirling Engine 

In [3] an acceptable dwell mechanism is expressed, able to be used in a a-Stirling 
engine. But this one is not good enough regarding the phases of the two pistons 
(Fig. 3). In order to improve this situation it is necessary to change the position of 
the articulation B between bars 2 and 3. It is preserved the relative position of the 
cylinder axis parallels. Hence the mechanism will become like the schema pre¬ 
sented in Fig. 4. 

The same reference xoy axes are placed in accordance with the Schmidt Analysis 
(Gustav Schmidt has made the classical Stirling Analysis in 1871). 

Using the same reference we can replace the pseudo-sin function of compression 
cylinder, in Schmidt analysis by the new space function, in order to develop a new 
isothermal analysis of the engine. 


3 On the Movement Functions of the Mechanism 

The purpose of the analysis is to find out the movement functions yC(t) and yE(t) to 
reveal if there is a proper phase relation between them or not. The proper phase 
relation demonstrates that the mechanism is good for the application. Second, it 
must show that yE(t) has a large enough stop level, in accordance with the 
expansion stroke time of the expansion cylinder. 


Fig. 3 General bar 
mechanism 
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Fig. 4 The modified bar mechanism 


The new displacement function yE(t) will be inserted in Schmidt analysis, 


instead of sin relation. 

Notations used: 

01A = Re rod 
AC = b; BB’ = h; AB’ = bl 
02D = R2 
BD = K 

The point A(xA, yA) has a position expressed by the equation system (1). The 
angle 0 = 0(t) and 0o = Oo when the rod 01A is oriented to oy axis. 


%a Rcrod sin(0) 
y A = h\ +R C rod cos (0) 


( 1 ) 


Using the sin theorem equation (2) can be written as follows: 


b 


Rc rod 
sin(a) 


( 2 ) 


sin(7i — 0) 


resulting the angle a = a(t): 



a = arc sin 


b 


( 3 ) 
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The point B(x b , y b ) has a position expressed by the equation system: 

j Xb = (b\ — b ) sin a — h cos a 
\ yb = hi + R c rod cos (0) — b\ cos a — h sin a 



where bi, b, h are the geometrical parameters of the crank (Fig. 5). 

The point C(x c , y c ) has a position expressed by the equation system: 

f X c = 0 

\ y c = h\ + R c rod cos(d) -b cos a 


The point D(x d , y D ) has a position expressed by the equation system: 


J X[) = e — R i sin /? 

\ yD = h 2 + Ri cos p 


The point E(x e , y E ) has a position expressed by the equation system: 


( x E = e 

\ yE = h 2 - 2 R 2 cos p 
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In order to solve y E (t) the expression of (3(t) is needed. To obtain the function 
(3(t) the equation system of the BD element must be written: 

J x B + K cos 5 + R 2 sin p — e = 0 
( y B — K sin 5 + R 2 cos P — h 2 = 0 


The equation system (8) has two variables: (3 and 5. The variable 5 is not needed, 
so it can be eliminated. 


j x B + R 2 sin P — e = —K cos 5 
\ ys + R 2 cos P — h 2 = K sin S 


2 and + 



resulting the equation: 

^ ^2 2R 2 ^ b ~ &) sin ft 

+ 2R 2 {y B ~ h 2 ) cos ft - 2x B e - 2y B h 2 = K 2 


( 10 ) 


It must be applied the next substitution: 

A = 2R 2 (x B - e) 

B = 2R 2 (y B — h 2 ) 

C = K 2 + 2x B e + 2y 5 /z 2 - 

2 d2 2 72 

-yB~ R 2~ e ~ h 2 

The equation (10) becomes simpler: 


A sin p + B cos p = D 

(ID 

A\/l — cos 2 P + B cos P = D 

(12) 

„ _ D 2 B 2 BD 

cos p rcos p - 2—^ cos p 

A 2 A 2 A 2 

(13) 

(B 2 \ BD D 2 

\A 2 + 1 2 A 2C ° SlS + A 2 1_0 

(14) 



The second substitution will be: 
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resulting an algebraic equation: 


PX 2 -QX + R = 0 


(15) 


where X = cos /3 

The equation that remains, gives the (3 values (16) because (3 = (3(0) or (3 = (3(t). 


/? = arc cos 


Q+^Q 2 +4PR 
2 P 


(16) 


The roots with negative sign are not practical for the mechanism; they are 
however in the field of arcos function. 

Replacing (3 into y E expression provides the function y E = y E (t). There is 
impossible to obtain y E = y E (t) as an analytical expression. Hence, we need to 
calculate y c and y E position point by point. 

Using a computer program the moving functions y c = y c (t) and y E = y E (t) can be 
obtained. The two functions are shown in Fig. 6, y c as y-expansion and y E as 
y-compression. 
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Fig. 6 Piston strokes 
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4 The Schmidt Analysis Insertion Paragraph 


From the Schmidt analysis only the paragraphs where it is necessary to replace the 
old values by the new ones are shown. The original calculus presented in [5] is 
preserved and marked. 

Therefore the mass of gas inside the regenerator space is given by: 

(56) m r = p.V r .ln(T h /T c )/(R(T h - T c )) 

The volume of the compression space varies in sinusoidal fashion: 

(57) V c = V clc + 0.5.V SWC (1 + cos(0)) 

0 = crank angle, proportional to time 

V c ic = clearance volume = minimum of V c 

Vswc = swept volume, max{ V c } = V c i c + V swc 

The volume of the expansion space varies in sinusoidal fashion: 

(58) V e = Vcie + 0.5.V swe (l + cos(0 + 5)) 


The equation (S8) has to be improved by replacing: O.5.V swe .cos(0 + 5) with y E (t). 
The same operation has to be done with the equation. (SI3) in Schmidt Analysis. 


(S13) p 


m to tRT h _1_ 

V re f (e + c + d) + c cos 6 + e cos (6 + 3) 


Equation (S13) allows to determine the pressure p at any angle 6 for given 
volumes and temperatures T h and T c . 

By differentiating equation (S8) with respect to 6 , the following equation from 
paragraph 3.1 in [5] are obtained: 


dVe = —O.5.V swe .sin(0 + 5).d0 


Expression 0.5.V swe .sin(6 + 3).d6 ” must be replaced with a numeric differ¬ 
entiation of yE(t). 

One has to note that the sense of the crank rotation is in accordance with the 
classical Schmidt analysis calculus. 


5 Conclusion 

The low efficiency of an a-Stirling engine drives us to improve on its work 
characteristics in order to enhance its efficiency. The main way to achieve this is 
to improve on the engine characteristics [6] . There are a lot of mechanisms species 
that could be used to obtain a stop movement necessary for the compression 
cylinder. A bar mechanism could be more reliable than a cam mechanism despite 
that the cam mechanism is more precise. 

The bar mechanism presented in this paper has seven bars but seems to be 
reliable enough for this application. 
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There are, also, other ways to improve the efficiency, like changing the working 
gas in the cylinders. A good efficiency is obtained using a gas that has a high 
specific heat, like hydrogen or helium. 
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